NASA  TECHNICAL  NOTE 


I 


D^STRfBUTIOM  STATEi^^lEMT  A 
Approved  for  Public  Release 
Distribution  Unlimited 

MONTE  CARLO  APPROACH 
TO  TOUCHDOWN  DYNAMICS 
FOR  SOFT  LUNAR  LANDING 

by  Robert  E.  Lavender 

George  C.  Marshall  Space  Flight  Center 
Huntsville,  Ala. 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION  -  WASHINGTON,  D.  C.  -  NOVEMBER  1965 


NASA  TN  D-3117 


MONTE  CARLO  APPROACH  TO  TOUCHDOWN 

DYNAMICS  FOR  SOFT  LUNAR  LANDING 

By  Robert  E.  Lavender 

George  C.  Marshall  Space  Flight  Center 
Huntsville,  Ala. 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


For  sale  by  the  Clearinghouse  for  Federal  Scientific  and  Technical  Information 
Springfield,  Virginia  22151  -  Price  $3.00 


TABLE  OF  CONTENTS 


Page 


SUMMARY .  1 

INTRODUCTION  .  1 

VEHICLE  CHARACTERISTICS .  1 

INITIAL  CONDITIONS .  4 

TOUCHDOWN  DYNAMICS  PROGRAM  .  9 

RESULTS .  11 

APPENDIX  A.  TOUCHDOWN  CONDITIONS .  37 


APPENDIX  B.  RANDOM  NUMBER  GENERATOR  CONFIDENCE  TEST  ....  47 


LIST  OF  ILLUSTRATIONS 

Figure  Title  Page 

1.  Vehicle  Configurations .  17 

2.  Spring  Constants  for  Vehicles  I,  E,  and  El  .....  j, .  18 

3.  Force-Stroke  Relationship .  19 

4.  Chi-square  Frequency  Function .  20 

5.  Chi-square  Distribution  Function . 21 

6.  Normal  Distribution  Function  .  22 

7.  Surface  and  Body  Axes  Orientation .  23 

8.  Symmetric  Heading  Orientations .  24 

9.  Stability  Boundary  for  Vehicle  3 .  25 

10.  Stability  Boundary  for  Vehicle  7 .  26 

11.  Lower  Limit  Probability  of  Stable  Landing . * .  27 

iii 


LIST  OF  TABLES 


Table  Title  Page 

1.  Geometric  Properties  of  LFV,  LEM,  and  LLS .  28 

2.  Geometric  Properties  of  Vehicles  I,  n,  and  El .  28 

3.  Limiting  Forces  for  Type  I  Vehicles .  29 

4.  Dynamic  Scaling  Factors .  29 

5.  Initial  Conditions  from  Chi-Square  Distribution .  30 

6.  Tabulation  of  Distribution  Functions .  31 

7.  Summary  of  Touchdown  Dynamics  Runs  for  Vehicle  1  32 

8.  Summary  of  Touchdown  Dynamics  Runs  for  Vehicle  5  34 

9.  Summary  of  Cases  which  Tumble .  35 

10.  Sample  Probability  of  Stable  Landing .  36 

11,  Lower  Limit  Probability  of  Stable  Landing  .  36 


IV 


DEFINITION  OF  SYMBOLS 


Symijol 

C 

CSA 

C 

nip 

C 

sp 

D 

O' 

fc> 

k 

V 

k 

z 

L 

111 

n^, 

11 

N 


r 


V 

V 


Definition 
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MONTE  CARLO  APPROACH  TO  TOUCHDOWN 
DYNAMICS  FOR  SOFT  LUNAR  LANDING 

SUMMARY 


^"^esults  of  analytical  touchdown  dynamics  investigations  are  presented 
which  were  conducted  to  obtain  estimates  of  the  probability  of  stable  landing  for 
configurations  with  various  landmg  gear  diameters.  Both  three-legged  and  four 
legged  vehicles  are  considered  in  the^analysis.  Dynamic  scaling  considerations 
are  taken  into  account  so  that  the  results  are  applicable  to  a  wide  range  of 
vehicle  size  and  mass.  A  Monte  Carlo  app^ach  is  taken  in  the  determination 
of  initial  landing  conditions,  j  /  ^ 

I  Results  indicate  that  for  a  given  probability  of  stable  landing,  a  three- 
legged  vehicle  requires  a  landing  gear  diameter  only  slightly  larger  than  the 
diameter  required  for  a  four-legged  vehicle.  Therefore  the  three-legged 
vehicle's  landing  gear  should  weigh  ^ 


INTRODUCTION 


In  the  overall  design  and  performance  of  spacecraft  intended  to  soft-land 
on  the  moon,  analysis  of  touchdown  dynamics  motion  is  an  important  part.  Such 
analysis  in  the  past  has  been  restricted  to  specific  missions  and  vehicles  designed 
for  those  missions  such  as  Surveyor  (ref.  1),  Lunar  Excursion  Module  (ref.  2), 
Lunar  Flying  Vehicle  (ref.  3),  and  Lunar  Logistics  System  (ref.  4).  In  general, 
the  analysis  for  each  vehicle  has  not  been  a  complete  parametric  study  over  a 
wide  range  of  initial  conditions,  but  rather  a  limited  study  with  emphasis  on  what 
the  vehicle  can  do  under  the  worst  set  of  conditions .  The  parametric  approach  is 
not  feasible  because  of  the  large  number  of  parameters  involved.  For  example, 
if  a  parametric  study  of  a  particular  vehicle  is  considered  with  only  three 
variations  of  ten  initial  touchdown  conditions,  a  total  of  three  to  the  tenth  power, 
or  59,  049  combinations,  is  produced. 

Presented  in  this  report  are  results  of  a  general  touchdown  dynamics 
study  in  which  a  Monte  Carlo  approach  is  taken  to  the  determination  of  initial 
conditions.  This  approach  is  more  realistic  than  simply  choosing  the  worst-case 
conditions  and  does  not  require  the  prohibitive  amount  of  analysis  which  a 


parameti'ic  approach  would  require.  The  vehicle  configurations  used  are 
determined  in  a  standardized  form,  and  the  results  are  equally  applicable  to  the 
Lunar  Flying  Vehicle,  Lunar  Excursion  Module,  or  Lunar  Logistics  System,  even 
though  these  vehicles  vary  greatly  in  size  and  mass. 

Initial  conditions  have  been  determined  using  a  Monte  Carlo  approach  for 
400  cases.  For  each  case,  the  initial  conditions  were  determined  for  10  parame¬ 
ters:  vertical  velocity,  horizontal  velocity,  velocity  cross-slope  angle,  vehicle 
pitch,  bank,  and  heading  angles,  lunar  slope,  and  the  vehicle *s  pitch,  roll,  and 
yaw  rates  at  touchdown.  The  coefficient  of  friction  between  the  landing  pads  and 
lunar  surface  is  assumed  to  be  infinite  (no  sliding) .  Eight  generalized  vehicles 
are  analyzed  for  these  400  cases  to  determine  whether  or  not  the  vehicle  being 
considered  lands  safely  or  tumbles.  The  eight  vehicles  considered  correspond 
to  three-legged  and  four-legged  vehicles  with  four  variations  in  the  ratio  of 
landing  gear  diameter  to  center-of-gravity  height.  From  these  results,  an 
estimate  is  made  of  the  probability  of  stable  landing  on  the  lunar  maria  for  each 
of  the  eight  vehicles. 

The  author  is  indebted  to  Mr.  John  D.  Capps,  Computation  Laboratory, 
who  programmed  the  equations  for  the  three-dimensional  touchdown  dynamics 
digital  pi’ogram  as  well  as  the  program  used  in  the  determination  of  initial 
conditions. 


VEHICLE  CHARACTERISTICS 


The  Lunar  Flying  Vehicle  (LFV)  is  a  small  vehicle  which  has  received 
consideration  for  lunar  surface  exploration  after  the  initial  Apollo  landings. 
Primary  emphasis  has  been  placed  upon  a  vehicle  having  an  80-kilometer 
maximum  range  when  carrying  a  payload  of  two  men  with  their  portable  life 
support  systems.  The  primary  mission  for  such  a  vehicle  was  considered  to  be 
rescue  back  to  the  Lunar  Excursion  Module  (LEM)  from  a  disabled  surface 
roving  vehicle.  A  secondary  mission  was  that  of  supplementing  the  roving 
vehicle  by  permitting  flights  into  areas  which  are  inaccessible  to  a  roving  vehicle. 

The  Lunar  Logistics  System  (LLS)  is  a  large  cargo-carrying  vehicle 
which  has  received  study  as  a  means  of  delivering  15-ton  payloads  to  the  moon  in 
support  of  a  lunar  base. 

Even  though  the  three  spacecx^aft  (LFV,  LEM,  LLS)  vary  greatly  in  size 
and  mass,  all  three  vehicles  can  be  reduced  to  approximately  the  same  nondimen- 
sional  geometry  as  shown  in  Table  1,  As  a  result  of  this,  three  standard  vehicles 
(I,  II,  and  m)  have  been  established  which  are  representative  of  the  LFV,  LEM, 
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and  LLS  size  vehicles.  The  standardized  dimensions  are  shown  in  Table  2.  The 
radius  of  gyration  about  the  bank  (roll)  axis  is  assumed  equal  to  that  shown  for 
pitch.  The  radius  of  gyration  in  yaw  (about  the  vehicle's  centerline)  is  taken  to 
be  0.9  times  the  pitch  value.  Having  standardized  the  three  vehicles,  all  touch¬ 
down  dynamics  calculations  performed  for  Vehicle  I  can  be  made  applicable 
through  dynamic  scaling  factors  to  Vehicle  n  and  Vehicle  HI. 

Eight  vehicles  are  considered,  as  shown  in  Figure  1,  corresponding  to 
four  values  of  the  ratio  of  landing  gear  diameter  to  center-of-gravity  height.  Both 
three-legged  and  four-legged  vehicles  are  considered  even  though  the  LFV,  LEM, 
and  LLS  have  all  received  major  emphasis  as  four-legged  vehicles.  A  major 
purpose  of  this  study,  however,  is  to  obtain  estimates  of  the  relative  probability 
of  stable  landing  with  three-legged  vehicles  compared  to  four-legged  vehicles. 
Although  the  diameter — c.g.  height  ratios  for  LFV,  LEM,  and  LLS  vary  from 
2.  50  to  3. 25  (Table  1) ,  the  values  chosen  for  this  generalized  study  vary  from 
2.0  to  3.0.  The  large  landing  gear  diameter  designed  for  the  LLS  was  based 
upon  very  severe  initial  conditions  which  are  believed  to  be  overly  conservative, 
resulting  in  an  excessively  heavy  landing  gear. 

The  spring  rates  assumed  for  the  vehicles  of  this  study  (Fig.  2)  are  the 
equivalent  vertical  spring  rates  per  leg  when  landing  on  a  level,  high-friction 
surface  with  no  sliding  of  the  foot  pads.  These  vertical  spring  rates  were 
obtained  by  assuming  an  elastic  spihng  rate  (7005  N/ cm  or  4000  lb/ in  for 
Vehicle  I)  along  the  three  struts  of  an  inverted  tripodal  leg.  As  seen,  the 
vertical  spring  rate  decreases  as  the  landing  gear  diameter  increases,  negating 
somewhat  the  beneficial  result  of  the  larger  landing  gear  diameters.  The  spring 
rates  assumed  for  the  LLS  size  vehicle  (Type  III)  are  rather  large  and  may  be 
very  difficult  to  obtain.  These  spring  rates  are  required  to  satisfy  dynamic 
similarity.  The  results,  of  this  study  may  therefore  be  somewhat  optimistic  for 
this  larger  vehicle  because  of  the  high  stiffness  assumed.  On  the  other  hand, 
the  assumed  spring  rate  for  Vehicle  I  (LFV)  is  probably  low  so  that  the  results 
obtained  may  be  somewhat  pessimistic.  Previous  study  (ref.  5)  has  shown  that 
decreasing  stiffness  has  a  detrimental  effect  on  landing  stability. 

The  force-stroke  relationship  in  each  strut  is  shown  in  Figure  3.  Upon 
compression,  each  strut  is  assumed  to  compress  elastically  until  the  limiting 
force  in  compression  is  reached.  Further  stroking  takes  place  with  a  constant 
stroking  force  assumed.  As  the  strut  begins  to  re-extend,  the  stored  elastic 
energy  is  released  until  the  strut  force  reaches  zero.  Further  re-extension  of 
the  strut  takes  place  under  zero  load  or  free-return  conditions.  No  limiting 
force  in  tension  is  assumed.  The  same  spring  constant  is  assumed  for  the  strut 
in  tension  as  used  in  compression.  The  limiting  force  is  assumed  to  represent 
the  crushing  force  of  aluminum  honeycomb  energy  absorbers  in  the  struts.  Upon 
subsequent  stroking,  the  strut  shortens  under  zero  load  as  long  as  it  is  still  in  the 
free-return  stroking  region. 
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The  limiting  force  or  crushing  force  values  for  the  eight  vehicles  were 
determined  so  that  the  deceleration  load  factor  is  limited  to  3.  0  earth  g’s  when 
the  vehicle  lands  on  a  level,  high  friction  surface  on  all  legs  simultaneously.  This 
load  factor  appears  to  be  a  good  choice  from  previous  study  (ref.  6) .  If  the  chosen 
load  factor  is  too  high,  a  landing  gear  designed  for  such  loads  is  too  heavy.  If  the 
load  factor  is  too  low,  excessive  stroking  is  required  to  absorb  the  landing  energy, 
thus  increasing  the  strut  length  required  and  landing  gear  weight.  Each  leg  is 
considered  to  be  an  inverted  tripodal  arrangement  of  one  main  strut  and  two  lower 
support  struts.  The  limiting  force  of  each  support  strut  is  considered  to  be  one- 
third  of  the  main  strut  limiting  force.  This  relationship  is  believed  to  be  a  good 
choice  based  upon  previous  analysis  (ref.  7) .  The  limiting  forces  obtained  for  the 
eight  vehicles  are  shown  in  Table  3.  These  values  correspond  to  the  forces  for 
Vehicle  I  (LEV) .  Corresponding  forces  for  the  larger  vehicles  are  increased  by 
the  dynamic  similarity  force  factors  given  in  Table  4. 

The  mass  of  Vehicle  I  (LEV)  is  considered  to  be  453.  59  kg  (earth  weight 
of  4448  N  or  1000  lb) .  The  mass  ratios  given  in  Table  4  can  be  used  to  obtain  the 
mass  of  the  larger  vehicles. 


INITIAL  CONDITIONS 


Initial  conditions  have  been  obtained  using  a  Monte  Carlo  approach  which 
randomly  selects  values  from  zero  to  one  and  relates  these  values  to  the  initial 
conditions  sought  through  the  distribution  functions  that  have  been  established  for 
the  initial  variables.  Determining  a  reasonable  distribution  function  for  each  of 
the  variables  of  interest  is  the  most  questionable  part  of  the  entire  study  because 
of  the  lack  of  data.  However,  one  value  of  this  study  may  be  that  of  helping  to 
point  out  the  need  for  detailed  and  realistic  statistical  information  concerning 
initial  touchdown  conditions. 

Various  simulations  have  been  performed  (refs.  8-13)  which  have  yielded 
some  information  on  the  conditions  at  touchdown  for  lunar  landing.  All  of  these 
simulations  have  been  conducted  with  varying  degrees  of  realism  relative  to 
instrumentation  and  control  response,  pilot  training,  test  objectives,  pilot  motiva¬ 
tion,  and  otlier  factors.  Results  must  be  viewed  with  some  degree  of  caution  since 
they  may  not  realistically  simulate  an  actual  vehicle  with  well  trained  and 
motivated  pilots.  Euture  efforts,  it  is  hoped,  will  provide  more  realistic  informa¬ 
tion  on  touchdown  conditions  through  the  use  of  the  Lunar  Landing  Research  Eacility 
(ref.  14)  and  also  the  free  flight  Lunar  Landing  Vehicle  (ref.  15) , 

As  far  as  the  vertical  and  horizontal  components  of  touchdown  velocity  are 
concerned,  it  is  reasonable  to  expect  that  the  vast  majority  of  cases  of  actual 
landings  with  well  motivated  pilots  will  be  very  soft  landings  with  low  touchdown 
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velocities.  The  simulations  reported  by  Hill  (ref.  9)  show  that  the  vertical 
descent  rate  is  less  than  1.2  m/s  (4  fps)  for  90  percent  of  the  cases.  The  forward 
velocity  at  touchdown  was  also  reported  as  less  than  0.6  m/ s  (2  fps)  for  90  per¬ 
cent  of  the  cases.  The  simulation  results  reported  by  Wood  and  Post  (ref.  12) 
indicate  more  severe  landing  velocities.  About  35  percent  of  the  time  the  pilots 
failed  to  land  within  the  design  envelope  of  3.  0  m/ s  (10  fps)  vertical  velocity, 
and  1.  5  m/s  (5  fps)  horizontal  velocity.  It  was  concluded  that  this  large  failure 
rate  can  be  attributed  not  only  to  difficulty  of  the  task,  but  also  to  the  training 
level  of  the  pilots  and  limitations  in  the  simulation.  It  was  believed  that  greater 
pilot  experience,  and  improvements  in  displays  and  display  precision  will  reduce 
such  failures  considerably  in  future  studies. 

It  is  expected  that  the  frequency  functions  of  the  vertical  and  horizontal 
velocity  components  should  have  some  tail-off  toward  higher  velocities  to  represent 
more  severe  cases.  The  frequency  function  chosen  as  a  reasonable  representation 
of  the  vertical  and  horizontal  velocity  components  is  given  in  Figure  4.  This 
function  corresponds  to  the  chi-square  frequency  function  with  three  degrees  of 
freedom.  Appropriate  choices  for  multiplication  constants  convert  the  values  of 
the  variable  x  to  the  vertical  and  horizontal  velocities. 

Concerning  the  lunar  slope  expected  in  maria  regions  of  the  moon,  it  is 
exiDected  that  the  vast  majority  of  landings  will  be  on  low  sloped  terrain  with  only 
a  small  percent  of  the  landings  on  rougher  terrain.  It  may  well  be  argued  that 
manned  landings  on  rough,  highly  sloped  terrain  will  never  take  place.  However, 
for  the  small  LFV  exploration  system,  some  landings  may  well  be  on  steeper 
slopes  because  of  the  desire  to  land  near  interesting  geological  and  topological 
features.  In  addition,  if  lunar  dust  somewhat  clouds  the  pilot's  vision  during  the 
final  phase  of  landing,  it  is  possible  that  one  or  more  pads  may  impact  on  a 
protuberance  or  into  a  depression  which  would  increase  the  effective  landing  slope. 
In  addition,  the  pad-surface  interaction  during  impact  may  well  serve  to  increase 
the  effective  landing  slope.  Experimental  evidence  (refs.  16-17)  indicates  that 
the  bearing  strength  of  the  soil  is  not  the  sole  criterion  for  penetration.  There  is 
a  tendency  for  the  downhill  pads  to  penetrate  the  soil  more  deeply  than  the  uphill 
pads  during  a  landing,  and  thus  increase  the  effective  landing  slope.  Considering 
these  factors,  it  is  assumed  that  the  frequency  function  shown  in  Figure  4  is  also 
a  reasonable  approximation  of  the  effective  slope  to  be  encountered  upon  landing. 
Mason  and  colleagues  (ref.  18)  have  established  a  similar  highly  skewed  frequency 
function  for  a  representative  lunar  surface  profile  from  20  100-kilometer  traverses. 
Their  work  was  performed  to  establish  a  baseline  for  determining  lunar  roving 
vehicle  requirements.  An  appropriate  choice  of  a  multiplication  factor  converts 
the  variable  x  in  Figure  4  into  the  effective  lunar  slope. 

The  distribution  (cumulative  frequency)  function,  shown  in  Figure  5,  is 
used  in  the  Monte  Carlo  approach  to  obtain  vertical  and  horizontal  velocity  compo¬ 
nents  and  the  effective  lunar  slope.  The  values  given  in  Table  5  have  been  obtained 
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from  this  distribution  function  by  use  of  appropriate  multiplication  constants.  The 
multiplication  constants  were  chosen  on  the  basis  of  the  values  emphasized  in 
Table  5.  For  example,  McCauley  indicates  (ref.  19)  that,  in  the  1-10  meter 
scale,  the  median  slope  of  the  maria  is  about  5  degrees.  Thus,  with  the  assumed 
distribution,  20  percent  of  the  maria  has  a  slope  of  about  2  degrees  or  less, 

80  percent  is  within  10  degrees,  and  1  percent  is  24  degrees  or  more. 

For  the  LEM  size  (TyiDe  II)  vehicle,  it  is  assumed  that  99.  865  percent  of 
the  landings  will  have  a  vertical  velocity  of  3.  048  m/ s  (10  fps)  or  less,  and  a 
horizontal  velocity  of  1.  524  m/ s  (5  fps)  or  less.  Thus,  with  the  assumed 
distribution,  90  percent  of  the  landings  have  a  vertical  velocity  of  1.  188  m/s 
(3.  90  fps)  or  less,  and  a  horizontal  velocity  of  0.  594  m/ s  (1.  95  fps)  or  less. 
While  the  values  were  based  upon  the  LEM  design  values  of  10  fps  and  5  fps 
(ref.  12),  the  90  percentile  values  also  correspond  quite  well  with  the  simulation 
results  of  Hill  (ref.  9).  Thus,  it  is  considered  that  the  LEM  size  vehicle  (TyiDe 
II)  will  land  50  percent  of  the  time  with  a  vertical  velocity  of  0.448  m/ s  (1.47  fps) 
or  less,  with  a  horizontal  velocity  of  0.  224  m/ s  (0.  735  fps)  or  less. 

For  the  results  to  be  equally  applicable  to  the  LFV  size  (Tyi^e  I)  vehicle 
and  the  LLS  size  (Type  III)  vehicle,  the  velocities  (also  shown  in  Table  5)  were 
adjusted  by  dynamic  scaling  considerations.  For  example,  it  is  assumed  that 
90  percent  of  the  landings  for  the  LLS  (Type  III)  vehicle  will  be  with  a  vertical 
velocity  of  1.638  m/ s  (5.374  fps)  or  less.  For  the  LFV  (Type  I)  vehicle,  it  is 
assumed  that  90  percent  of  the  landings  will  occur  with  a  vertical  velocity  of 
0. 706  m/ s  (2.  32  fps)  or  less. 

The  dynamic  scaling  factors  given  in  Table  4  show  that,  for  dynamic 
similarity,  the  velocity  is  scaled  as  the  square  root  of  the  vehicle's  linear  dimen¬ 
sion  when  landing  in  the  same  gravity  field.  Thus,  for  larger  vehicles,  larger 
landing  velocities  can  be  tolerated.  While  this  indicates  that  the  larger  vehicles 
may  be  easier  to  land,  this  table  also  shows  that  the  spring  rate  for  dynamic 
similarity  is  proportional  to  the  mass  and  inversely  proportional  to  the  linear 
dimension  when  landing  in  the  same  gravity  field.  Thus,  since  the  mass  increases 
to  some  power  of  the  linear  dimension,  higher  spring  rates  are  required  for  the 
larger  vehicles  for  dynamic  similarity. 

The  horizontal  velocity  is  assumed  to  be  controlled  as  far  as  possible  to  be 
zero.  Therefore,  it  is  considered  that  the  direction  of  the  horizontal  velocity 
vector  at  touchdown  with  respect  to  the  principal  slope  direction  has  an  equal 
probability  of  being  in  any  direction.  Further  study  of  this  problem  may  show  an 
advantage  to  having  a  horizontal  velocity  of  sufficient  magnitude  to  guarantee  an 
uphill  landing  rather  than  a  smaller  residual  component  that  may  be  in  the 
downhill  direction.  The  assumption  of  trying  to  land  at  zero  horizontal  velocity, 
however,  should  certainly  be  applicable  to  unmanned  vehicles  such  as  the  LLS  or 
an  unmanned  LEM  used  as  a  LEM  truck.  It  is  also  assumed  that  the  vehicle’s 
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geometric  heading  orientation  with  respect  to  the  principal  slope  direction  has  an 
equal  probability  of  being  in  any  direction.  This  assumption  is  applicable  to 
unmanned  vehicles  and  also  to  manned  vehicles  where  vision  impairment  due  to 
lunar  dust  might  make  the  determination  of  the  direction  of  principal  slope 
difficult.  Under  clear  vision,  the  pilot  would  probably  orientate  the  vehicle’s 
heading  to  align  in  a  1-2-1  orientation  to  the  slope  for  improved  stability. 

The  pitch  and  bank  angles  at  touchdown,  referenced  to  the  horizontal,  are 
assumed  to  follow  a  normal  distribution  with  zero  mean.  Likewise,  the  vehicle's 
angular  rates  in  pitch,  roll,  and  yaw  are  also  assumed  to  follow  the  same  distribu¬ 
tion.  The  normal  distribution  with  zero  mean  and  standard  deviation  of  one  is 
shown  in  Figure  6.  Appropriate  multiplication  constants  convert  the  x  variable 
to  angles  and  angular  rates .  The  standard  deviation  for  pitch  and  bank  angles  is 
considered  to  be  2  degrees  and  is  applicable  for  Vehicles  I,  II,  and  III.  The 
standard  deviation  of  pitch,  roll,  and  yaw  rates  is  taken  to  be  1.60  degrees  per 
second  for  the  Type  n  (LEM  size)  vehicle.  Past  simulation  studies  show  that 
these  are  reasonable  assumptions.  For  dynamic  similarity,  the  cori’esponding 
standard  deviation  of  angular  rate  is  taken  to  be  2.  70  degrees  per  second  for  the 
Type  I  (LFV)  vehicle  and  1. 16  degrees  per  second  for  the  Type  III  (LLS) 
vehicle . 

The  10  variables  of  vertical  velocity,  horizontal  velocity,  lunar  slope, 
cross-slope  angle,  heading  angle,  bank  angle,  pitch  angle,  pitch  rate,  roll  rate, 
and  yaw  rate  are  assumed  to  be  uncorrelated. 

The  Monte  Carlo  procedure  for  obtaining  the  initial  conditions  can  be 
summarized  in  the  following  manner.  A  set  of  values  for  400  cases  is  obtained 
where  each  case  consists  of  obtaining  10  pseudo-random  numbers  from  a 
previously  developed  digital  random  number  generator  subroutine  for  a  rectangular 
frequency  function  (ref.  20) .  These  numbers  are  equal  to  or  greater  than  zero 
and  equal  to  or  less  than  one.  These  10  random  numbers,  RN^  through  RNjo,  are 
then  used  as  F  (x^)  through  F  (xjq)  to  obtain  the  values  for  x^  through  x^q  from 
Table  6. 

This  table  was  obtained  from  values  tabulated  by  Abramowitz  and  Stegun 
(ref.  21) .  Linear  interpolation  is  used  to  obtain  the  values  of  x^  to  Xjq-  The 
initial  conditions  are  then  computed: 


V,  =  0.  0605  x, 
ht 

(i) 

(2) 

=  0.  121  X2 

(3) 

V  =  V^g  k  '  V 

V  »  Z  V 

(4) 
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y  =  2.  113  X3 

(5) 

eSA  =  360  X4 

(6) 

V,  =  cos  eSA 

h  ht 

(7) 

V  =  -  V,  ,  sin  eSA 
z  ht 

(8) 

0  0  ==  360  X5 

(9) 

i/jfl  =  2  xg 

(10) 

0  =  2  X7 

P 

(11) 

^0  ^ 

(12) 

CO  ^2,7  Xg 

X  ® 

(13) 

CO  =  2.  7  Xq 

y  ^ 

(14) 

CO  —  2 ,  7  X  j  A 
z 

(15) 

The  values  were  obtained  corresponding  to  the  Type  I  (LFV)  vehicle  for 
convenience  in  performing  touchdown  dynamics  computations  since  the  time 
required  to  obtain  the  touchdown  motion  is  less.  A  computing  interval  of  0.  002 
sec  was  used  in  performing  the  touchdown  dynamics  computations.  The  same 
results  could  have  been  obtained  with  the  larger  vehicles  in  the.  same  computation 
time  by  increasing  the  computing  interval  by  the  appropriate  time  factors  given  in 
Table  4.  Tabulation  of  the  initial  conditions  for  the  400  cases  is  given  in 
Appendix  A. 

A  I’andom  number  confidence  test  was  performed  to  judge  the  adequacy  of 
the  pseudo-random  numbers  generated  by  the  digital  subroutine.  The  results  are 
given  in  Appendix  B.  Four  thousand  random  numbers  were  generated  corresponding 
to  400  cases  with  10  random  numbers  in  each  case,  RNj  through  RN^o,  used  to 
obtain  the  10  initial  conditions.  Therefore,  400  values  were  obtained  for  each  of 
the  10  random  numbers  RNj  through  RNk,.  The  range  from  zero  to  one  was 
divided  into  10  intervals  so  that  these  400  values  would  have  an  expected  frequency 
of  40  in  each  interval.  The  observed  frequencies  were  obtained  and  the  chi-square 
value  obtained.  These  chi-square  values  should  follow  a  chi-square  distribution 
with  9  degrees  of  freedom.  Therefore,  the  value  of  chi-square  could  be  expected 
to  exceed  12.2  (ref.  21)  for  20  percent  of  the  cases,  or  for  2  of  the  10  random 
numbers.  None  of  the  random  numbers  have  a  chi-square  value  which  exceeds 
the  12.  2  value.  Therefore,  the  values  generated  by  the  random  number  subroutine 
are  adequate. 
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TOUCHDOWN  DYNAMICS  PROGRAM 


The  touchdown  dynamics  program  used  for  this  study  corresponds  to  a 
three-dimensional  method  programmed  for  digital  computer  operations.  The 
method  was  developed  as  a  simplified  version  of  the  landing  dynamics  method  of 
analysis  for  the  LEM  described  by  Mantus,  Lerner,  and  Elkins  (ref.  22),  This 
simplified  version  is  restricted  to  landings  on  a  plane  sloping  surface  without 
protuberances  or  depressions  and  without  sliding  of  the  landing  pads  along  the 
surface.  The  surface  and  body  axes  orientation  used  is  shown  in  Figure  7.  The 
velocity  cross-slope  angle  is  the  angle  from  the  principal  slope  direction  to  the 
total  horizontal  velocity  vector.  The  vehicle's  heading  angle  defines  the  orientation 
of  the  vehicle's  landing  pads.  Zero  heading  angle  corresponds  to  a  1-2-1  pad 
orientation  for  a  four-legged  vehicle  and  to  a  1-2  pad  orientation  for  a  three-legged 
vehicle.  Heading  angles  of  ±45  degrees  correspond  to  a  2-2  pad  orientation  for 
the  four-legged  vehicle.  Heading  angles  of  ±60  degrees  correspond  to  a  2-1  pad 
orientation  for  the  three-legged  vehicle.  These  various  orientations  are  shown  in 
Figure  8. 

As  previously  stated,  the  initial  heading  orientation  has  been  taken  to  have 
equal  probability  of  being  in  any  direction  (360  degrees) .  However,  as  far  as  the 
touchdown  dynamics  program  is  concerned,  0,  90,  180,  270,  and  360  degrees 
heading  orientation  is  still  a  1-2-1  orientation  for  a  four-legged  vehicle. 

Similarly,  0,  120,  240,  and  360  degrees  heading  orientation  is  still  a  1-2  orienta¬ 
tion  for  a  three-legged  vehicle.  Consequently,  the  initial  heading  angle  obtained 
from  the  Monte  Carlo  analysis  has  been  converted  to  an  equivalent  heading  within 
±45  degrees  for  the  four-legged  vehicle  and  within  ±60  degrees  for  the  three-legged 
vehicle  for  input  into  the  touchdown  dynamics  program.  Thus: 


4>  Qi~  ^  0 
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(16) 

<^04  =  00- 
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(17) 
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(19) 

P  0i~  P  0  ~ 

360, 

315  <  00  -  360 

(20) 
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^  03  ^  0  “  ^"^^5 


180  <  00  -  300 


(23) 


003=^00-360,  300  <00  ^  360  (24) 

These  equivalent  heading  angles  for  the  400  cases  are  also  tabulated  in  Appendix  A. 

Each  run  begins  with  the  vehicle's  lowest  pad  just  touching  the  surface. 
Strokes  and  forces  along  the  struts  are  computed  with  subsequent  determination 
of  forces  and  moments  referred  to  the  body  axes.  The  moments  are  used  to 
obtain  the  angular  accelerations  about  the  body  axes  and  are  integrated  numeri¬ 
cally  to  obtain  the  angular  velocities  about  the  body  axes.  The  rate  of  change  of 
the  angles  0,  and  0  are  then  computed  from  expressions  containing  the 
angular  velocities  about  the  body  axes.  The  angular  accelerations  of  0,  0,  0  are 
then  obtained  by  equations  derived  from  differentiating  the  equations  for  0,  <1,  6. 
The  angles  <9,  0,  and  0  are  obtained  by  numerical  integration  using  the  computed 
rates  and  angular  accelerations  for  these  angles.  The  forces  computed  are 
converted  to  forces  along  the  inertial  (surface)  axes  and  the  vehicle's  transla¬ 
tional  accelei’ations  obtained.  These  accelerations  are  integrated  numerically 
to  obtain  the  velocity  and  position  of  the  vehicle's  center  of  gravity.  The 
computations  are  repeated  at  each  time  step  until  sufficient  time  has  elapsed  to 
determine  whether  or  not  the  vehicle  tumbles  or  is  stable. 

As  long  as  a  pad  is  off  the  surface,  it  travels  along  with  the  vehicle.  Upon 
reaching  the  surface,  the  pad  remains  at  its  impact  position  until  the  forces  on  the 
pad  normal  to  the  surface  become  negative.  The  negative  force  is  set  to  zero  and 
the  pad  then  lifts  from  the  surface  and  again  travels  along  with  the  vehicle  until  it 
imioacts  the  surface  again. 

Examples  of  stability  boundaries  which  have  been  obtained  by  this  program 
are  shown  in  Figures  9  and  10,  Figure  9  gives  the  stability  boundary  in  terms  of 
the  lunar  slope  versus  the  vehicle’s  heading  angle  for  a  three-legged  vehicle. 

Figure  10  gives  corresponding  results  for  a  four-legged  vehicle.  As  seen  from 
Figure  9,  stable  landings  can  be  made  on  lunar  slopes  from  6  degrees  to  18 
degrees  depending  upon  the  vehicle’s  heading  orientation  with  respect  to  the 
principal  slope  direction.  For  the  four-legged  vehicle  (Fig.  10),  stable  landings 
can  be  made  on  slopes  from  10  degrees  to  about  18  degrees.  Minimum  stability 
does  not  correspond  to  45  degrees  heading  (2-2  orientation)  but  rather  to  a  heading 
angle  of  about  30  degrees  for  this  case.  The  impact  velocity  values  shown 
correspond  to  the  small  (Tyi3e  I)  vehicle. 
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RESULTS 


With  the  three-legged  vehicle  having  the  smallest  landing  gear  (Vehicle  1), 
the  400  cases  were  reviewed  to  determine  by  inspection  which  cases  were  more 
likely  to  tumble.  Most  of  the  cases  are  obviously  stable  because  of  the  small 
touchdown  velocities,  combined  with  small  surface  slopes.  Three  of  the  cases 
were  considered  to  obviously  tumble  because  of  the  surface  slope  being  greater 
than  24  degrees.  Touchdown  dynamics  motion  was  then  computed  for  45  cases 
which  were  believed  to  have  a  substantial  probability  of  tumbling.  Of  these,  14 
cases  did  result  in  tumbling  motion. 

Based  upon  these  results,  touchdown  motion  was  computed  for  an  addi¬ 
tional  34  cases  which  were  believed  to  be  stable,  but  with  a  substantial  degree  of 
uncertainty.  Of  these,  three  cases  tumbled.  Finally,  an  additional  nine  cases 
were  chosen  which  were  thought  definitely  would  be  stable  but  lacked  the  high 
degree  of  certainty  desired.  Of  these,  none  tumbled.  Based  upon  these  results, 
a  final  review  was  made  of  the  remaining  cases  with  the  conclusion  that  all  would 
certainly  be  stable.  Therefore,  for  Vehicle  1,  20  cases  of  the  400  correspond  to 
tumbling  cases. 

A  summary  of  the  cases  used  for  the  touchdown  dynamics  computations 
for  Vehicle  1  is  given  in  Table  7.  It  is  observed  that  tumbling  resulted  for 
Case  43  with  a  lunar  slope  of  only  6.  9  degrees  while  Case  332,  with  a  lunar  slope 
of  19.2  degrees,  was  stable.  While  the  lunar  slope  is  a  very  important  parameter 
to  consider  in  touchdown  dynamics  analysis,  other  initial  conditions  are  also 
obviously  important.  Case  43  corresponds  to  a  near  1-2  landing  orientation  with 
a  relatively  high  vertical  impact  velocity,  while  Case  332  corresponds  to  a  near 
2-1  orientation  with  a  small  impact  velocity. 

Touchdown  dynamics  motion  was  then  computed  using  Vehicle  2  for  the  20 
cases  which  were  previously  determined  as  tumbling  cases.  All  stable  cases  for 
Vehicle  1  would  also  be  stable  cases  for  Vehicle  2  since  nothing  is  changed  except 
the  landing  gear  diameter,  which  is  larger.  Of  these  20  cases,  11  cases  tumbled 
including  the  3  cases  considered  obvious  for  Vehicle  1. 

Touchdown  motion  was  tlien  computed  using  Vehicle  3  for  the  11  cases 
which  tumbled  when  using  Vehicle  2.  Of  these,  four  cases  tumbled. 

Touchdown  motion  was  tlien  computed  using  Vehicle  4  for  the  four  cases 
which  tumbled  when  using  Vehicle  3.  Of  these,  none  tumbled. 

A  similar  analysis  was  performed  for  the  four-legged  vehicles.  Vehicle  5 
was  used  and  touchdowm  dynamics  was  computed  for  25  cases  which  were  thought 
to  have  a  substantial  probability  of  tumbling.  Of  these,  11  cases  tumbled. 
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Computations  were  then  performed  for  32  additional  cases  which  were  believed  to 
be  stable,  but  with  a  substantial  degree  of  uncertainty.  Of  these,  three  cases 
tumbled.  Finally,  an  additional  14  cases  were  considered  which  were  thought 
definitely  would  be  stable  but  lacked  the  high  degree  of  certainty  desired.  Of 
these,  none  tumbled.  Three  cases  (with  a  slope  of  over  24  degrees)  were 
considered  to  be  obviously  tumbling  cases.  Therefore,  for  Vehicle  5,  17  of 
the  400  cases  correspond  to  tumbling  cases. 

A  summary  of  the  cases  used  for  touchdown  dynamics  calculations  for 
Vehicle  5  is  given  in  Table  8.  No  case  tumbled  corresponding  to  a  lunar  slope  of 
less  than  15.  5  degrees,  while  no  case  was  stable  with  a  slope  greater  than  19. 3 
degrees.  This  is  a  slope  spread  of  only  3.  8  degrees  compared  to  a  slope  spread 
of  12.3  degrees  for  the  corresponding  three-legged  vehicle  (Vehicle  1).  While  7 
cases  tumbled  for  Vehiele  1  witli  a  lunar  slope  of  less  than  15.  5  degrees,  only  13 
cases  tumbled  with  lunar  slopes  of  15.  5  degrees  or  more  compared  to  the  17 
cases  for  Vehicle  5.  These  results  indicate  that  a  relatively  well  defined  critical 
slope  can  be  established  for  four-legged  vehieles  compared  to  three-legged 
vehicles. 

The  vehicle  with  the  next  largest  diameter  (Vehicle  6)  was  then  eonsidered. 
Touchdown  motion  was  computed  for  the  17  eases  with  the  result  that  7  cases 
tumbled  ineluding  the  3  cases  considered  obvious  for  Vehicle  5.  These  seven 
cases  were  then  eonsidered  using  Vehicle  7  with  the  result  that  three  cases 
tumbled.  Finally,  these  three  cases  were  considered  using  Vehiele  8  with  the 
result  that  none  tumbled.  The  cases  corresponding  to  tumbling  motion  are 
summarized  in  Table  9. 

A  summary  of  the  sample  probabilities  for  stable  landing  is  given  in 
Table  10.  The  sample  probability  of  stable  landing  varies  from  0.  950  for  the 
three-legged  vehicle  with  the  smallest  landing  gear  to  1. 000  for  the  vehicles  with 
the  largest  landing  gear.  The  lower  limit  probability  for  stable  landing  with  0.95 
and  0. 995  confidence  coefficients  for  the  binomial  population  has  been  taken  from 
published  tables  (ref.  23)  based  upon  the  sample  size  of  400.  These  results  are 
given  in  Table  11.  The  table  also  shows  the  lower  limit  probability  for  stable 
landing  based  upon  an  approximate  equation  given  by  Dalton  (ref.  24) .  This 
equation  is 


P. 


(l-C)/e 


[(3  +  20/4] 


/  (n-n^/2), 


(25) 


where 


d  =  0  for  n^  =  0 
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6  =  1  for  0. 

This  equation  is  seen  to  give  a  very  good  approximation  of  the  results  for  the 
binomial  population. 

The  lower  limit  probability  fox  stable  landing  for  the  binomial  population 
with  a  confidence  coefficient  of  0.  995  is  shown  in  Figure  11.  Since  neither 
Vehicle  4  nor  Vehicle  8  had  any  cases  which  tumbled,  touchdown  dynamics  runs 
were  also  made  for  an  additional  three-legged  vehicle  and  an  additional  four¬ 
legged  vehicle  with  a  diameter-c.  g.  height  ratio  of  2.  75.  The  three-legged 
vehicle  tumbled  for  Case  317.  This  case  has  the  highest  lunar  slope  (28.2 
degrees)  of  the  400-case  sample  and  the  vehicle  lands  in  a  near  1-2  orientation. 
The  resulting  probability  for  stable  landing  for  this  vehicle  is  also  shown  in 
Figure  11.  The  four-legged  vehicle  tumbled  for  Case  129  and  also  for  Case  317. 
Case  129  has  the  second  highest  lunar  slope  (27.  3  degrees)  of  the  400-case 
sample.  The  resulting  probability  for  stable  landing  is  also  shown  in  Figure  11. 

The  probability  values  for  the  three-legged  vehicles  vary  in  a  smooth 
manner,  whereas  the  values  for  the  four-legged  vehicles  are  somewhat  erratic. 
At  a  diameter-c.  g.  height  ratio  of  2.75,  the  probability  of  stable  landing  for  the 
four-legged  vehicle  is  less  than  for  the  three-legged  vehicle  based  upon  the  400- 
case  sample.  It  is  believed  that  this  result  would  not  be  obtained  if  the  sample 
size  were  increased  by  an  order  of  magnitude.  Similarly,  at  a  diameter-c.  g. 
height  ratio  of  2.25,  the  probability  of  stable  landing  for  the  four-legged  vehicle 
appears  to  be  somewhat  high.  Consequently,  the  curve  shown  in  Figure  11  for 
the  four-legged  vehicles  has  been  faired  in  a  manner  to  remove  the  erratic 
nature  of  the  data  points. 

The  results  indicate  that,  for  a  given  probability  of  stable  landing,  a 
three-legged  vehicle  requires  a  landing  gear  diameter  only  slightly  larger  than 
the  diameter  required  for  a  four-legged  vehicle.  Since  each  leg  is  folded  for 
flight  to  the  moon  and  must  be  deployed  before  landing,  the  reliability  of  landing 
gear  deployment  should  be  greater  for  a  three-legged  vehicle.  With  only  three 
landing  pads  to  contact  the  surface,  there  is  less  chance  of  contacting  a 
protuberance  or  depression.  A  three-legged  vehicle  assures  a  more  positive 
final  resting  support,  whereas  a  four-legged  vehicle  may  tend  to  rock  back  and 
forth.  With  only  a  slightly  larger  diameter  required,  the  three-legged  vehicle's 
landing  gear  should  weigh  less. 


George  C.  Marshall  Space  Flight  Center, 

National  Aeronautics  and  Space  Administration, 
Huntsville,  Alabama,  September  30,  1965 
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FIGURE  2.  SPRING  CONSTANTS  FOR  VEHICLES  I,  II,  AND  III 
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FIGURE  3.  FORCE-STROKE  RELATIONSHIP 
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FIGURE  4.  CHI-SQUARE  FREQUENCY  FUNCTION 
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FIGURE  5.  CHI-SQUARE  DISTRIBUTION  FUNCTION 
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FIGURE  7.  SURFACE  AND  BODY  AXES  ORIENTATION 
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CSA  =  30 ,  My-l.S  m/sec,  V^t  =  0.9  m/sec 
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FIGURE  10.  STABILITY  BOUNDARY  FOR  VEHICLE  7 


FIGURE  11.  LOWER  LIMIT  PROBABILITY  OF  STABLE  LANDING 
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TABLE  1.  GEOMETRIC  PROPERTIES  OF  LFV,  LEM,  AND  LLS 


V  ehicle 

D 

L 

Ya 

Y 

P 

r 

kz 

cm  (in) 

cm  (in) 

cm  (in) 

cm  (in) 

cm  (in) 

cm  (in) 

LFV 

279  (110) 

112  (44) 

74  (29) 

30  (12) 

56  (22) 

58  (23) 

LEM 

813  (320) 

320  (126) 

170  (67) 

119  (47) 

208  (82) 

1 65  (65) 

LLS 

1966  (774) 

605  (238) 

381  (150) 

180  (71) 

330  (130) 

259  (102) 

V  ehicle 

D/L 

Yp/L 

r  /L 

k  /L 
z 

LFV 

2.  50 

0.  66 

0.  27 

0.  50 

0.  52 

LEM 

2.  54 

0.  53 

0.  37 

0.  65 

0.  52 

LLS 

3.  25 

0.  63 

0.  30 

0.  55 

0.  43 

TABLE  2.  GEOMETRIC  PROPERTIES  OF  VEHICLES  I,  II,  AND  III 


V  ehicle 

Type 

Y; 

a/J- 

Y 

p/-*-' 

r 

/L 

I 

LFV 

0. 

61 

0, 

,  31 

0. 

57 

0. 

48 

II 

LEM 

0. 

61 

0, 

,  31 

0. 

57 

0. 

48 

III 

LLS 

0. 

61 

0, 

.  31 

0. 

57 

0. 

48 

Vehicle 

Type 

Y. 

a 

Y 

P 

k 

cm 

(in) 

cm 

(in) 

cm 

(in) 

cm 

(in) 

I 

LFV 

68 

(27) 

35 

(14) 

64 

(25) 

54 

(21) 

II 

LEM 

19^ 

i  (77) 

99 

(39) 

182 

>  (72) 

153 

(60) 

III 

LLS 

369  (145) 

188 

(74) 

345  (136) 

290 

(114) 
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TABLE  3.  LIMITING  FORCES  FOR  TYPE  I  VEHICLES 


V  ehicle 

N 

D/L 

Cmp 

c 

sp 

ky 

N 

lb 

N 

lb 

N/  cm 

lb/  in 

1 

3 

2.  00 

3576 

804 

1192 

268 

8756 

5000 

2 

3 

2.  25 

3897 

876 

1299 

292 

7425 

4240 

3 

3 

2.  50 

4218 

948 

1406 

316 

6252 

3570 

4 

3 

3.  00 

4977 

1119 

1659 

373 

4588 

2620 

5 

4 

2.  00 

2682 

603 

894 

201 

8756 

5000 

6 

4 

2.  25 

2922 

657 

974 

219 

7425 

4240 

7 

4 

2.  50 

3162 

711 

1054 

237 

6252 

3570 

8 

4 

3.  00 

3738 

840 

1246 

280 

4588 

2620 

TABLE  4 

.  DYNAMIC  SCALING  FACTORS 

Vehicle 

Type 

Quantity 

Scale  Factor 

U 

III 

Acceleration 

g/gj  =  CL 

1. 00 

1. 00 

Mass 

m/m^  ^  p 

14.  00 

42.  00 

Liength 

.  L/Lj  =  y 

2.86 

5.41 

Angular  Acceleration 

cl/ y 

0.  35 

0.  18 

Time 

(y/ci) 

1.  69 

2.  33 

Angular  Velocity 

(a/y)’^' 

0.  59 

0.43 

Angle 

(a/y)“ 

1. 00 

1. 00 

Velocity 

(a  .  y) 

1. 69 

2.  33 

Force 

a  .  p 

14.00 

42.00 

Mass  Density 

P/y^ 

0.  60 

0.  27 

Moment  of  Inertia 

P  •  y" 

115  1229 

Stress 

a  •  p/y^ 

1.71 

1. 44 

Spring  Rate 

n  *  p/y 

4.  90 

7.  76 
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TABLE  5.  INITIAL  CONDITIONS  FROM  CHI-SQUARE  DISTRIBUTION 


F  (x) 

% 

ZO. 

000 

50. 

000 

80. 

000 

X 

- 

1. 

005 

Z. 

366 

4. 

642 

\t 

- 

0. 

061 

0. 

143 

0. 

281 

v 

0. 

12Z 

0. 

286 

0. 

562 

V 

y 

deg 

Z. 

I 

5. 

0"= 

9. 

8 

Vht 

I 

> 

m/  s 

0. 

057 

0. 

133 

0. 

261 

Vht 

,II 

m  /  s 

0. 

096 

0. 

224 

0. 

442 

\t 

,  III 

m/  s 

0. 

133 

0. 

310 

0. 

608 

Vv, 

I 

m/  s 

0. 

114 

0. 

266 

0. 

522 

II 

m/  s 

0. 

192 

0. 

448 

0. 

884 

III 

m/  s 

0. 

266 

0. 

620 

1 . 

216 

90. 

000 

95. 

000 

99. 

000 

99. 

865 

99. 

997 

6. 

251 

7. 

815 

1 1 . 

341 

16. 

022 

Z1 . 

000 

0. 

378 

0. 

473 

0. 

686 

0. 

970 

1. 

271 

0. 

756 

0. 

946 

1. 

372 

1 . 

939 

Z. 

541 

13. 

2 

16. 

5 

Z4. 

0 

33. 

9 

44. 

4 

0. 

353 

0. 

441 

0. 

637 

0. 

901 

1. 

181 

0. 

594 

0. 

744 

1 . 

078 

1. 

524-:= 

1 . 

996 

0. 

819 

1 . 

0Z4 

1 . 

486 

Z. 

100 

Z. 

752 

0. 

706 

0. 

88Z 

1. 

274 

1 . 

802 

z. 

362 

1. 

188 

1 . 

488 

Z. 

1  56 

3. 

048=:= 

3. 

992 

1. 

638 

Z. 

048 

z. 

972 

4. 

200 

5. 

504 

=!=V;’lu,es  used  to  establish  multiplication  constants. 
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TABLE  6.  TABULATION  OF  DISTRIBUTION  FUNCTIONS 


Normal  Distribution  Chi-Square  Distribution 


o 

1! 

cr  =  1) 

(V  =  3) 

F(X) 

X 

F  (X) 

X 

0. 00000 

-5.  00 

0 

0 

0. 00003 

-4.  00 

0. 08111 

0.  50 

0. 00024 

-3.  50 

0.  19875 

1. 00 

0. 00135 

-  3.  00 

0. 31773 

1.  50 

0. 00621 

-2.  50 

0. 42759 

2.  00 

0. 02275 

-  2.  00 

0. 52444 

2.  50 

0. 04006 

-1.75 

0. 60837 

3.  00 

0. 06681 

-1.  50 

0. 67900 

3.  50 

0.  1 0565 

-  1.  25 

0.  73854 

4.  00 

0.  1  5866 

-  1. 00 

0. 82820 

5.  00 

0. 22663 

-  0.  75 

0. 88839 

6.  00 

0.  30854 

1 

o 

o 

0. 92810 

7.  00 

0. 40129 

-  0.  25 

0. 95399 

8.  00 

0. 50000 

0 

0. 97071 

9.  00 

0. 59871 

0.  25 

0. 98143 

10.  00 

0. 69146 

0.  50 

0. 99262 

12.  00 

0. 77337 

0.  75 

0. 99818 

15.  00 

0. 84134 

1.00 

0. 99956 

18.  00 

0. 89435 

1.  25 

0. 99997 

21.  00 

0. 93319 

1.  50 

1 . 00000 

O 

o 

ro 

0. 95994 

1.75 

0. 97725 

2.  00 

0. 99379 

2.  50 

0. 99865 

3.  00 

0. 99976 

3.  50 

0.  99997 

4.  00 

1. 00000 

5.  00 
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TABLE  7.  SUMMARY  OF  TOUCHDOWN  DYNAMICS  RUNS  FOR  VEHICLE  1 


Group  A 

Group  B 

Group  C 

Case 

Result 

y,  deg. 

Case 

Result 

y.  deg. 

Case 

Result 

y,  deg 

6 

S 

13.  0 

145 

s 

10.  4 

13 

s 

14.  8 

8 

S 

10.  3 

168 

s 

12.  7 

61 

s 

9.  5 

9 

S 

11.8 

172 

s 

1.  7 

95 

s 

00 

14 

s 

14.  2 

235 

s 

10.  0 

100 

s 

14.  3 

39 

s 

15.7 

238 

s 

9.7 

102 

s 

13.  9 

43 

T 

6.  9 

242 

s 

12.  4 

103 

s 

4.9 

46 

s 

00 

247 

s 

00 

123 

s 

00 

• 

o 

64 

s 

12.  4 

250 

s 

o 

00 

243 

s 

12,2 

66 

s 

1,2.  8 

259 

s 

6.  8 

337 

s 

8.  3 

82 

s 

15.  3 

275 

s 

13.  4 

84 

T 

15.  5 

277 

s 

15.  0 

87 

T 

11.5 

288 

s 

12.  4 

118 

T 

10.  7 

289 

s 

00 

00 

120 

S 

10.  3 

291 

s 

7.2 

144 

S 

14.  1 

296 

s 

8.  6 

148 

s 

11.8 

303 

s 

9.  5 

157 

T 

23,  0 

307 

s 

14.8 

1  58 

S 

16.  0 

309 

s 

10.7 

173 

T 

15.  5 

314 

s 

13.  0 

177 

T 

19.  3 

315 

s 

12.  1 

179 

S 

10.  4 

319 

s 

7.  3 

187 

s 

9.  1 

321 

s 

13.  5 

198 

s 

14.  9 

327 

T 

14.  9 

206 

s 

16.8 

3  38 

s 

6.  9 

207 

s 

16,  9 

341 

s 

8,  4 

21 1 

s 

11,3 

357 

s 

15.  3 

212 

T 

9.7 

363 

T 

14.  0 

214 

S 

11.  3 

364 

s 

15.  4 
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TABLE  7  (CONTINUED) 


Group  A 

Group  B 

Case 

Result 

y>  ^eg. 

Case 

Result 

y>  deg. 

215 

s 

17.7 

368 

s 

14.  1 

234 

s 

14.  4 

369 

T 

12.  5 

294 

s 

16.  4 

386 

S 

9.  5 

299 

s 

11.6 

396 

S 

13.  3 

328 

s 

17.7 

397 

S 

8.  7 

332 

s 

19.2 

400 

S 

6.7 

333 

s 

16.  3 

339 

T 

17.  0 

345 

S 

17.2 

346 

T 

21.8 

Group  A: 

Cases  which  were  thought  to  have  a 

350 

T 

20.  6 

substantial  probability  of  tumbling. 

351 

T 

16.  5 

Group  B: 

Cases  which  were  believed,  with 

371 

T 

21.2 

substantial  uncertainty,  would  be 

379 

S 

14.  7 

stable. 

389 

S 

10.  9 

Group  C: 

Cases  which  were  thought,  with  only  a 

390 

S 

18.  3 

small  degree  of  uncertainty,  would 

398 

T 

23.  5 

be  stable. 
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TABLE  8.  SUMMARY  OF  TOUCHDOWN  DYNAMICS  RUNS  FOR  VEHICLE  5 


Group  A 

Group  B 

Group  C 

Case 

Result 

y,  deg. 

Case 

Result 

deg. 

.  Case 

Result  y,  deg. 

14 

S 

14.  2 

6 

S 

13.  0 

13 

S 

14.8 

39 

S 

15.7 

9 

S 

11.8 

64 

S 

12.  4 

66 

S 

12.  8 

43 

S 

6.9 

87 

S 

11.5 

84 

S 

15.  5 

46 

S 

8.9 

100 

S 

14.  3 

157 

T 

23.  0 

61 

S 

9.  5 

102 

S 

13.9 

173 

T 

15.  5 

82 

S 

15.  3 

118 

S 

10.7 

177 

S 

19.3 

95 

S 

8.  4 

145 

S 

10.  4 

198 

S 

14.  9 

103 

S 

4.9 

148 

S 

11.8 

206 

T 

16.8 

123 

S 

8.  0 

275 

S 

13.4 

215 

S 

17.7 

144 

S 

14.  1 

303 

S 

9.  5 

289 

S 

8.8 

158 

S 

16.  0 

321 

S 

13.  5 

296 

S 

8.  6 

207 

T 

16.  9 

338 

S 

6.9 

315 

S 

12.  1 

212 

S 

9.7 

389 

S 

10.  9 

319 

s 

7.  3 

228 

S 

9.  3 

396 

S 

13.  3 

328 

T 

17.7 

234 

s 

14.  4 

333 

S 

16.  3 

238 

s 

9.7 

339 

T 

17.0 

270 

s 

6.  7 

346 

21.8 

277 

s 

15.  0 

350 

T 

20.  6 

288 

s 

12.  4 

351 

T 

16.  5 

291 

s 

7.  2 

364 

S 

15.  4 

294 

T 

16.  4 

368 

S 

14.  1 

307 

S 

14.  8 

371 

T 

21.2 

314 

S 

13.  0 

Group  A: 

Cases  which  were 
thought  to  have  a 

390 

T 

18.  3 

327 

S 

14.  9 

substantial  proba¬ 

398 

T 

23.  5 

332 

T 

19.2 

bility  of  tumbling 

337 

S 

8.  3 

Group  B: 

Cases  which 

were 

345 

S 

17.2 

believed,  with  sub¬ 
stantial  uncertainty, 

357 

S 

15.  3 

would  be  stable 

363 

S 

14.  0 

Group  C: 

Cases  which 

were 

369 

S 

12.  5 

thought,  with  only 

376 

S 

5.  4 

a  small  degr 
uncertainty, 

ee  of 

would 

379 

S 

14.  7 

be  stable 
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TABLE  10.  SAMPLE  PROBABILITY  OF  STABLE  LANDING 


Ps 

:::  (400  - 

n^)  /400 

Vehicle 

1 

2 

3 

4  5 

6  7  8 

Ps 

0. 9500 

0.  9725 

0. 9900 

1.0000  0.9575 

0.9825  0.9925  1.0000 

TABLE  11.  LOWER  LIMIT  PROBABILITY  OF  STABLE  LANDING 


Vehicle 

C  -  0.  950 

1  2 

3 

4 

5 

6 

7 

8 

p  (Ref  Z3) 

0.  9282 

0.  9549 

0. 9773 

0.  9925 

0. 9369 

0.  9674 

0.  9807 

0.  9925 

p  (Ref  24) 
s 

0. 9303 

0.  9575 

0. 9787 

0.  9925 

0. 9394 

0. 9696 

0.  9817 

0.  9925 

V  ehicle 

C  =  0.  995 

1  2 

3 

4 

5 

6 

7 

8 

Pg  (Ref  23) 

0. 9149 

0. 9439 

0. 9689 
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APPENDIX  B 

RANDOM  NUMBER  GENERATOR  CONFIDENCE  TEST 


(RANDOM  NUmBEP  QEMtHATOF?  CONFIDENCE  TPST  TOR  A  FLAT  D I STR I  BUT  I  biM  . 
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HAHDOM  number 

“TTpTFT 

•ATUR  CONFIDENCE  TEST 

FOlP 

A  flat oTSTRlbUtlU^. 
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number  n^Vt-ATUP  COMFjDPNCE  TFST  FOP  A  FLAT  D  j  STR  i  bUT  j  , 
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RAMDOM  number  HFM'fcHATUW  CONFIDFNCF;  TPST  Fn^  A  PLAT  D I  SIR  I  bUT  1  U^' . 
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RANDOM  number  RFMfcf-ATUW  COMFjDFNnF  TPST  FDR  4  FLAT  D  j  STR  j  bUT  j  ON  . 
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